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Abstract 

A novel method of summation for power series is developed. The method is based 
on the self-similar approximation theory. The trick employed is in transforming, first, a 
series expansion into a product expansion and in applying the self-similar renormaliza- 
tion to the latter rather to the former. This results in self-similar factor approximants 
extrapolating the sought functions from the region of asymptotically small variables 
to their whole domains. The method of constructing crossover formulas, interpolat- 
ing between small and large values of variables is also analysed. The techniques are 
illustrated on different series which are typical of problems in statistical mechanics, 
condensed-matter physics, and, generally, in many-body theory. 
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1 Introduction 



It is the standard situation for practically all realistic problems in physics and applied math- 
ematics that their solutions are obtained by some kind of perturbation theory or iterative 
procedure resulting in asymptotic series in powers of a parameter or a variable. Such series 
are usually only asymptotic and diverge for finite values of the expansion variable, while what 
is of physical interest are exactly the finite values of the latter. That is why the problem of 
defining an effective sum of divergent asymptotic series is of paramount importance. The 
most often employed summation techniques are Pade and Borel summations [1,2]. These 
techniques have a number of limitations, because of which in many cases they are either of 
bad accuracy or often are not applicable at all. 

An alternative approach to the problem of defining effective sums of asymptotic series or 
effective limits of iterative sequences has been developed, being based on the ideas of optimal 
control theory, renormalization group theory, and general dynamical theory. The principal 
concepts of this approach are as follows. 

The first pivotal step is the introduction of control functions whose role is to govern an 
optimal convergence of approximation sequences [3]. This results in the optimized pertur- 
bation theory which has been widely employed for a variety of applications (see surveys in 
[4,5]). 

The second idea is to consider the passage from one successive approximation to another 
as the motion on the manifold of approximants, where the approximant order plays the role 
of discrete time. The recurrent relations, representing this motion, are formalized by means 
of group self- similarity [6-11]. 

A dynamical system in discrete time, whose trajectory is bijective to an approximation 
sequence, is called the approximation cascade. Embedding the latter into an approximation 
flow makes it possible to derive differential and integral equations of motion, whose fixed 
points represent self-similar approximants to the sought function [6-11]. The stability of the 
calculational procedure is characterized by local multipliers [8-11]. 

Another key point of the approach is the introduction of control functions with the help 
of fractal transforms [12-14]. This allows one to analyse asymptotic series by transforming 
them into their fractal counterparts and then invoking all the machinery of the self-similar 
approximation theory for the transformed series. In this way, the approximants, possessing 
a nice self-similar structure, were obtained, such as the self-similar exponential approximants 
[13-15] and self-similar root approximants [16-18]. 

In the present paper, we suggest and give a detailed analysis of one more type of approx- 
imants possessing a nice self-similar structure, which, because of their form, we name the 
self- similar factor approximants. In Section 2, we explain the origin of these approximants 
and discuss some variants of their usage. Then, in Section 3, we illustrate the method on 
several functions containing exponentials. Section 4 contains examples of more complicated 
and more realistic physical problems. Quartic anharmonic models are treated in section 5. 
Calculations of critical points and critical indices for some critical phenomena are presented 
in Section 6. The results are summarized in Section 7, where further possible developments 
of the method are discussed. 
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2 Method of self-similar factor approximants 



Suppose we are looking for a physical quantity described by a real function f{x): M — >^ M. It 
is the common case for realistic physical problems that the sought function is defined through 
so complicated equations that the problem can be solved only approximately, by invoking a 
kind of perturbation theory. Assume that the latter yields perturbative approximants 

fix) ^ Mx) (x^O), (1) 

of order A: = 0, 1,2,..., in the asymptotic vicinity of a; = 0. For a while, the concrete meaning 
of the function f{x) and of the variable x is not important. In particular applications, x may 
be a spatial variable, a coupling parameter, or any other physical parameter. This may also 
be a time variable in expansions having to do with differential evolution equations [19-21]. 

We may assume, without loss of generality, that the sought function f{x) has the property 
/(O) = 1, so that /fc(0) = 1. As is evident, this can be achieved for any function by the 
appropriate scaling and normalization. Then an asymptotic series in powers of x takes the 
form 

k 

/fc(a;) = l + $]a„xV (2) 

n=l 

If we apply the technique of the self-similar approximation theory to the fractal transform 
of series (2), wc would come, depending on the choice of control functions, to the self-similar 
exponential or root approximants [13-18]. But here, we follow another route. 

We may notice [22] that each series (2) can always be equivalently presented as a product 

k 

Mx) = f[{l + bix) , (3) 

i=l 

with the coefficients bi defined by the accuracy-through-order relations with respect to series 
(2). Now, we apply the procedure of self-similar renormalization not to the whole series (2) 
but separately to each factor of the product (3). This results in the transformation 

l + bix^{l + Aix)""' , (4) 

in which Ai and rii are control parameters (see details in [12-18]). In this way, starting with 
product (3), we come to the self-similar factor approximants 

f:{x)^l[{i+Axr . (5) 

i=l 

Control parameters Ai and rii can be found by means of the rc-cxpansion procedure [23], 
expanding expression (5) in powers of x and equating the like terms from series (2). It 
is also possible to define some of the controllers by matching the behaviour of the factor 
approximant (5) with the large-x asymptotic behaviour of f{x), provided this is available. 
For instance, if the asymptotic form 

fix) ~ /oo [x ^ oo) (6) 
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is known, then we may require the vahdity of the crossover condition 

f[Ar^foo, jZnr^P. (7) 

i=i 1=1 

The resulting crossover formula (5) would sew together two asymptotic expressions, at small 

X — > and large x oo. 

As a zero-order approximant, we may set 

f*(x) = l + aix. (8) 

The first-order approximant (5) is 

fl{x)^{l + A^xr, (9) 

and for the second-order, we have 

/*(x) = (l + Ar(l + A2xr . (10) 

Thus, we construct approximants of any order. For each approximant fl{x), the accuracy- 
through-order relation defines its own controllers Ai — Aik and rii — riik- However, to avoid 

cumbersome notation, we shall not employ the double indexation, keeping in mind that the 
controllers Ai and rii are a sort of variables that are different for different approximation 
orders k = 1,2, . . .. 

Each self-similar factor approximant (5) of order k contains 2k controllers Ai,...,Ak 
and ni, 712, rife. To define these through the re-expansion procedure, one needs to have 
2k nontrivial terms of series (2). Thus, an even-order perturbative series (2) uniquely defines 
the corresponding factor approximant (5). The possibility of using odd-order series will be 
discussed in Section 7. 

It may happen that, calculating the controllers A^ and rii, we shall obtain for some of 
them complex values. If the sought function is real, then the complex-valued controllers 
should come in complex conjugate pairs, so that approximant (5) be real. This means that 
the latter may contain the factors of the type with both z and a being complex valued. 
Since 

_ |-;|RcQ exp { — (Im a) arg 

is a nonalgebraic function of z, hence expression (5) may, in general, be a nonalgebraic 
function. According to definition [24], an analytical function that is not algebraic is termed 
transcendental. 

In this way, formula (5) can represent the following functions. When all powers rii are real 
integers, positive or negative, then fk{x) is a rational function. In the case when all powers 
are equal to ±1, expression (5) coincides with the form of Pade approximants, which are, 
therefore, just a narrow particular case of the self-similar factor approximants. When some 
of rii are real but not integer, then /^(x) is an irrational function. And if some of Ai and rii 
are complex valued, then /^(x) represents a transcendental function. As far as any sought 
function, having the same form (5), is exactly reproducible by /^(x), whose controllers are 
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defined by the re-expansion procedure [23], this means that the class of functions, that can 
be exactly reconstructed by the self-similar factor approximants, is rather wide, including 
rational, irrational, and some transcendental functions. This class is essentially wider than 
that related to Pade approximants and consisting solely of rational functions. Consequently, 
the factor approximants (5) can provide much better accuracy than Pade approximants for 
the majority of functions and also can be employed when the latter are not applicable at all. 

Approximants (5) can well approximate such entire transcendental functions as exponen- 
tials. This can be easily understood noticing that if a factor (l + Ax)^ is such that A = a/n, 
where a is a constant and n — ^ oo, then such a factor gives 



that is, it is equivalent to an exponential. Therefore, if in some particular calculations we 
find a very small A but a very large n, this hints that the corresponding factor (1 + Ax)"' 
is close to the exponential exp{nAx). This is why the factor approximants can provide a 
reasonable approximation for exponential functions. 

3 Reconstruction of functions including exponentials 

Exponential functions are ubiquitous is natural sciences representing various physical laws. 
At the same time, this is a simple example of an entire transcendental function. The best 
way of approximating these functions would be by means of the self-similar exponentials 
[13-15], which often reconstruct exponential- type functions exactly [25]. But the factor 
approximants, as is explained in the previous section, should also mimic well exponential 
behaviour. This is illustrated in the following examples. 

3.1 Nonmonotonic exponential function with maximum 

We choose a nonmonotonic function in order to stress that such a nonmonotonic behaviour 
presents no problem for the method. Let us consider the function 



Its asymptotic expansion at a; — > yields the finite sums (2) with the coefficients 
Qi = 1 , as = -2 , ag = 2.042 , 04 = -0.271 , 

05 = -3.572 , ag = 8.636 , ay = -12.661 , Og = 12.183 , 

and so on. To approximate function (11), we construct the self-similar factor approximants 
(5). The first-order approximant (9) is yet too simple to catch well the peculiarities of 
the function variation. So, here and in what follows, we shall start the analysis with the 
second-order approximant (10). 




{n — > 00) , 




(11) 
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Here for f2{x), we have 

Ai = 1.658 + 0.840 i, A2 = Al ; 

ni = -0.244- 1.077 i , 712 = . 

But, as a; — > 00, the function /2 (a^) decreases to zero, since ni+n2 — —0.488, while function 
(11) tends to unity. 
For we get 

^1 = 1.3540, ^2 = 1.0015 + 0.4034 i , A^ ^ A^ ; 

ni = 11.0051 , n2 = -5.4336 + 3.7402 i, 713 = 71^. 

Note that the enumeration of the controllers Ai and Ui can be arbitrary, as far as the order 
of factors in Eq. (5) is interchangeable. As a rule, we shall enumerate them according to the 
descending order of the values |Re Ai\. Now, rii + n2 + = 0.138, and the approximation 
is not bad till x ~ 150. 
For we find 

^1 = 1.1455, ^ = 1.0076 + 0.1703 i , ^3 = , ^4 = 0.5156 ; 

711 = 84.0147, 7i2 = -40.1038 + 36.6118 i , 713 = 71^ , 714 = -3.7834 . 

The sum of the powers decreases, 7ii + 712 + 713 + 714 = 0.024. And the accuracy is very good 
up to X of order of hundred, with the error being less than 5%. 



3.2 Nonmonotonic exponential function with minimum 

Let us now turn to the function 

/W-^{-(T^) (12) 

having a minimum. The corresponding coefficients in series (2) are 

Qi = -1 , 02 = 3 , 03 = -7.042 , 04 = 15.354 , 

05 = -32.261 , 0,6 = 65.951 , ar = -131.714 , ag = 257.734 . 

This function also tends to 1 at x — > 00. 
For we have 

Al = 1.65790 + 0.83984 i, A2^Al; 

ni = 0.24381 + 1.07665 i , n2 = n\. 

This approximant works well till x ~ 10, but at x ^ 00 its behaviour is not correct since 
ni+n2 = 0.488. 
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The controllers of are 

yli = 1.35398, = 1.00149 + 0.40337 i , A3 = ; 

ni = -11.00514 , 712 = 5.43355 - 3.74022 i , m^n\. 

This function decreases at x — > 00, since ni + ^2 + 713 = —0.133. Reasonable accuracy is 
provided till x ~ 10. 

The next approximant fXk^^ the parameters 

Ai = 1.14550 , A2 = 1.00759 + 0.17032 i , A3 = A^ , A4 = 0.51563 ; 

ni = -84.01466 , 712 = 40.10379 - 36.61184 i , 723 = , 714 = 3.78340 . 

The function decreases much slower than the previous approximations, since 7Zi + 772 + 713 + 
774 = —0.024. The accuracy is good up to x ~ 100. For instance, the percentage error at 
X = 40 equals -2%. 

3.3 Spectral density of black-body radiation 

This is given by the Planck formula p{x) — Cx^jie^ ~ !)• This is also a nonmonotonic 
function with a maximum. To get the presentation (2), we define j{x) = p{x)/Cx'^, which 
results in 

/(^) = ^ , (13) 
— 1 

where x = hw/kBT. The coefficients of series (2) are 

1 1 n ^ 

Ol = - 7: , ^2 = T7: , 03 = , 04 = 



2 ' 12 ' " ' 720 ' 

1 1 

05 = , ae = , 07 = , as = - 



30240 ' ' ' » 1209600 
_ _ 1 _ 691 

Og — , Oio — , Oil — , ai2 



47900160 ' - ' - 1307674368000 

For /I {x) , we get 

^1 = -0.05 + 0.119024 i , A2 = Al- 

771 = -1 + 2.520504 i , 772 = 77*. 

The function decreases at x — > 00, since 77i + 772 = —2, but not so fast as Eq. (13). 
For /Kx), one has 

= -0.052201 , A2 = -9.613987 x 10"^ + 0.150710 7 , A3 = A; ; 

m = 7.654768 , 773 = -1.327382 + 0.417823 7 , 773 = 77^ . 

The function increases at x — > 00 because 77i + 772 + 773 = 5, which is not correct. 
The controllers of fl{x) are 

Ai = -0.026326+0.044214 i , A2 = A* , A3 = -1.452250x10-^+0.158034 i , A4 = A* ; 
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ni = -0.939416 + 5.998214 i , ni^n\, rig = -1.060584 + 0.070005 i , 714 = . 

The sum ni + 712 + ^3 + ^4 = —4 shows that the function decreases as x — > 00, in agreement 
with the behaviour of Eq. (13). 

To analyse the behaviour of higher-order approximants, we have also calculated f^{x), 
for which J2l=ink = 0.367295; because of this f^{x) increases at infinity, though not as fast 
as fsix). And for feix), we find I]fc=i^fc ~ —3.300722; hence this approximant decreases 
as X — > 00. With increasing order k, the interval of x, where fkix) provides a very good 
approximation to f{x), becomes larger. Figure 1 presents the difference /K^) — f{x). 

3.4 Specific heat of diatomic gas 

The oscillational part of the specific heat for a diatomic gas has the form [26] 

Ci-) - j^^, , (14) 
with X = fiuj/kBT. The expansion coefficients of series (2) are 



Oi — — — , a-2 — —— , as — — , 04 — 



12 ' ' 240 ' ' 6048 ' " 172800 
1 691 



^5 — ~ r^r,r,r,Ar, 1 



5322240 ' ^ 118879488000 ' 
and so on. Expression (14) tends to zero as a; ^ oo. 
For the second approximant C2{x), we find 

Ai = 0.0251 , A2 = 0.0026 ; 

ni = -2.0656 , 712 = -11.9344 . 

The approach to zero, as x — > 00, is sufficiently fast, since 7^l + 7^2 = —14. 
For C3 (x) , one gets 

Ai = 0.025148 , A2 = 198.278139 , A3 = 2.629939 x 10 

ni = -2.065650 , na = , = -11.934396 . 

Here, ■0,1+77,2+^3 = —14. Both Cl.{x) are of good accuracy, with the difference C^{x) — C{x) 
less than 10^"^. Note that small A3 with large 123 hint at the exponential dependence. The 
comparison of 6*3(0;) with C{x), given in Eq. (14), shows that this factor approximant 
is essentially more accurate than the best Fade approximants P[4/2](a;) and P[5/i](a;), as is 
illustrated in Fig. 2. 



-3 
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3.5 Regular solution from Debye-Huckel theory 

Consider the function 

fix) - 1 (l- e-) (15) 

appearing in the Debye-Huckel theory [26]. This function contains an exponential, but the 
large-x asymptotic behaviour is of power law, 

fix) ~ 2 f- - ^) ix^oc). 



X x^ 



The coefficients of series (2) are 



1 



^^-"3' ^^-12' ^^-"60' "'- 360' 
11 11 



^5 — ~ oror^' ^6 — 7^7777:7: ) 07 — — — -— -TT) ^8 — 



2520' " 20160 ' 181440' 1814400 

Constructing the factor approximants (5), we shall invoke the crossover condition (7). 
For f2ix), we have 

Ai = 0.177337 + 0.176694 i , A^^ A\\ 

m = -0.5 + 0.441432 i , 712 = < . 

The maximal error of this approximant is about 2.5% at x~ 15. 
The controllers of f^ix) are 

Ai = 0.153026 , A2 = 0.076787 + 0.156796 i, A^^A^; 

m = -1.146332 , 712 = 0.073166 + 0.539398 i , nz = n\. 

The maximal error now is around —1% at x ~ 20. Percentage errors of f^ix) and the 
two-point Pade approximant P[3/2](a;) are shown in Fig. 3. 
For jl ix) , we get 

Ai = 0.095338 + 0.056647 i , ^2 = ^1 , A3 = 0.043027 + 0.139489 i , ^4 = ^3; 
m = -0.966783 + 0.967772 i , ris = n* , n^ = 0.466783 + 0.285025 i , ^4 = n* . 

Here the maximal, with respect to all x G [0, cxd), error is about 0.5% at a; ~ 30. The best 
two-point Pade approximant [1] involving the same number of terms has the maximal error 
—6% at X 20, which is an order larger than the maximal error of fXi^)- The accuracy of 
the latter may be improved even more, if to take into account two terms from the large-x 
expansion. Then the maximal error can be reduced to 0.15%. 



4 Structure factors and radiation intensity 

The quantities such as structure factors and radiation intensity are directly observable in 
experiments. That is why much attention is paid to their theoretical description. Below we 
show how one can derive rather simple and quite accurate expressions for such quantities, 
by means of the self-similar factor approximants. 
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4.1 Structure factor of branched polymers 

The structure factor S{x) is usually expressed in scaled units as a function oi x = k'^B? /d, 
where /c is a wave vector, R is gyration radius, and d is spatial dimensionality [27,28]. For 
(i = 3, the structure factor of branched polymers is given [28] by a confluent hypergeometric 
function S{x) = Fi(l; 3/2; 3x/2). The behaviour at large x is 

I 

S{x) ~ - x^^ {x — > oo) . 

At small X, the coefficients of series (2) are 

Qi — —1 , a2 = 0.6 , as = —0.257 . 
We shall use the crossover condition (7). Then for (x), we find 

^1 = 0.247387 + 0.483295 i , A2 = Al ; 

m = -0.5 - 0.778627 i, ria = < . 
The approximant S2{x) agrees well with the exact S{x), as is shown in Fig. 4. 

4.2 Structure factor of ring polymers 

The structure factor of three-dimensional dilute ring polymers S{x) can be considered as a 
function oi x = {kRY, where k is wave vector and R is gyration radius [29]. In series (2), 
we have 

ai = - - , 02 = 0.0610 , 03 = -0.0073 . 

Knowing the large-x asymptote 

S{x) ~ 0.63 x-^/^'' {u = 0.588) , 

we shall employ the crossover condition (7). Then we get for S2{x) 

yli = 0.065483 + 0.185055 i , A2 = Al ; 

ni = -0.425170 + 0.750184 i, 7x2^ n\. 
This gives a rather good accuracy, as compared to numerical results [29]. 

4.3 Luminescent intensity of donor-acceptor recombination 

The expression for the temporal dependence of the luminescence intensity, arising from 
radiative recombination of donor-acceptor pairs, has the form [30,31] 

7(z) = -exp{aX(z)}-^X(^) , 
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in which z = wt, w is a, transition probabihty per unit time, t is time, a is a constant, and 



K{z) = [exp (-^e-^) - 1 

We shall consider the reduced quantity 



dx 



with the typical parameter a — 6.514 x 10"'''. Then the coefficients of short-time series (2) 
are 

ai = -0.125 , 02 = 0.019 , ag = -2.604 x 10"^ , 
a4 = 3.334 X 10"^ , 05 = -3.858 x 10"^ , ae = 4.050 x 10"^ . 
In the long-time limit, one has 



/(^)-^^ exp|--(ln^)3| 



[z ^ 00] 



which tends to zero. 
For /2*(-z) we find 



Ai = 0.132653 + 0.072661 i, A2 = Al ; 



ni = -0.256845 + 0.391265 i , n^^nX. 

This approximant tends to zero, as 2; — > 00, since \ n2 — —0.514. 
For jl{z\ we obtain 

^1 = 0.098681 , A2 = 0.087209 + 0.090440 i, A3 = A* ; 

ni = -1.292027 , ria = 0.293510 + 0.269217 i , ^ n* . 

This also tends to zero, as 2; — > 00, with the power law rii + n2 + — —0.705. In the 
range < ^ < 10^, the maximal error is about 2% at 2; ~ 200. The best Fade approximant 
P[3/3]{z) is much less accurate, as is seen in Fig. 5. 



5 Models of quartic anharmonic oscillators 

Anharmonic oscillator models, especially low- dimensional ones, are usually considered as 
touchstones for checking the validity of any approximation technique. This is because these 
models mimic the main mathematical features of more realistic, but much more complicated, 
problems in many-body theory and field theory. Below, we illustrate the applicability of the 
method of self-similar factor approximants to two low-dimensional anharmonic models. 



11 



5.1 Partition function of zero-dimensional model 

The partition function is modelled by the integral 

Z{g) = — exp (-(^^ - gcp^) dip , (16) 

in which g imitates a coupling parameter. Coefficients of series (2) are given by the formula 

(-1)" ^ / 1\ 

In the strong-coupling limit, one has 

Z{g) ~ 1.022765 g-^l^ (g ^ oo) . (17) 
Following the technique of Section 2, we find for ^l(g') the controllers 

Ai = 19.141 , A2 = 4.859 ; 
m = -0.00862 , 712 = -0.120 . 

At large g, this yields 

Z2*(^7)^ 0.806 ^-^^-^^^ (^^00). 

The approximant gives a reasonable accuracy at finite g and is qualitatively correct at large 
g, though the decrease of the partition function (16) is not as fast as in Eq. (17). The 
amplitude of the decay is predicted with an error —21% and the power, with an error 48%. 
For Zl{g), we get 

Ai = 31.220 , A2 = 13.317 , A3 = 3.464 ; 

m = -0.000526 , 712 = -0.022 , 713 = -0.125 . 
In the strong-coupling limit, 

Zl((7)=. 0.807 {g^oo). 

Here the prediction errors are —21% for the amplitude and 41% for the power. 
The approximant Z^{g) contains 

Ai = 43.965 , A2 = 23.064 , A3 = 10.294 , A4 = 2.677 ; 

m = -0.000025 , 7i2 = -0.00259 , 713 = -0.035 , 714 = -0.124 . 
The strong-coupling limit gives 

Z:(^)~ 0.810 ^-"^-^^^ (^^00). 

The errors are —21% for the amplitude and 36% for the power index. 



12 



For Zl{g), we find 

Ai = 57.165211 , A2 = 33.720143 , A3 = 18.543443 , 

Ai = 8.396565 , A5 = 2.174639 ; 
ni = -1.025651 x 10"^ , ria = -2.150967 x 10"^ , rig = -6.028095 x 10"^ , 

714 = -0.044164 , 715 = -0.119596 . 
The strong- couphng hmit is 

0.81445 (5^00). 

Here, the errors of the amphtude is —20% and that of the power index is 32%. 
And for Zl{g\ we obtain 

Ai = 66.087061 , A2 = 41.445220 , A3 = 25.096672 , 

A4 = 13.752049 , A5 = 6.067701 , Aq = 1.511098; 
7ii = -1.153412 X 10"^ , 7i2 = -3.352991 x 10"^ , 713 = -1.330855 x 10~^ , 
7i4 = -0.014434 , 7i5 = -0.057625 , Uq = -0.110557 . 
In the strong-couphng hmit, this leads to 

0.82548 ^-°-i84 (^^00). 

The errors are —19% for the amphtude and 26% for the index. These results demonstrate 
monotonia numerical convergence. 

The accuracy of approximants can be essentially improved by invoking the crossover 
condition (7). Then for ^2(9')) we have 

Ai = 8.57593 , A2 = 0.30987 ; rii = -0.08136 , 712 = -0.16864 ; 

A^^ = 1.023 , m + 7i2 = -0.25 . 

This approximant provides a good accuracy in the whole region oi g & [0, 00) , with the 
maximal error of 6% at ~ 5. 
For Z^{g), we find 

Ai = 19.47265 , A2 = 5.21916 , A3 = 0.18018 ; 

ni = -7.94804 x 10"^ , 712 = -0.10947 , 713 = -0.13258 ; 

A^^ A^' = 1.023 , m + 7i2 + na = -0.25 . 

The maximal error is about 4% at ~ 10. 
The next approximant Z1[g) contains 

Ai = 31.45590 , A2 = 13.55695 , A3 = 3.72340 , A4 = 0.12535 ; 
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ni = -4.96209 x 10"^ , n2 = -0.02119 , rig = -0.11631 , 714 = -0.11200 ; 

A^^ A^^ A2^ = 1.023 , ni + 712 + 713 + 714 = -0.25 . 

The maximal error is around 3.5% a,t g ^ 10. 
And for Z^{g), we obtain 

Ai = 44.148724 , A2 = 23.247563 , A3 = 10.482932 , 

Ai = 2.878801 , A5 = 0.095481; 

7ii = -2.395229 x 10"^ , 713 = -2.477078 x 10"^ , 713 = -0.033250 , 

7i4 = -0.115813 , 7i5 = -0.098436 . 

Now the maximal error is close to 3% at gf ~ 20. 
Finally, for Zl{g), we have 

Ai = 57.315168 , A2 = 33.870044 , A3 = 18.694023 , 

^4 = 8.552264 , A^ = 2.338677 , Aq = 0.076825 ; 
m = -9.879942 x 10"^ , na = -2.072431 x 10"^ , 7x3 = -5.811092 x 10"^ , 
n4 = -0.042533 , 715 = -0.112765 , 7Z6 = -0.088682 . 
For this approximant, the maximal error is about 2.5% at ~ 30. 

5.2 Ground-state energy of the anharmonic oscillator 

We consider the one-dimensional anharmonic quartic oscillator with the Hamiltonian 

where x e (— cxd,+cxo). Perturbation theory gives [32] the ground-state energy E{g) as the 
series 

00 

^(^) = E«n^", (19) 
n=0 

where the first eleven coefficients are 

oo = 0.5 , oi = 0.75 , 02 = -2.625 , 03 = 20.8125 , 04 = -241.2890625 , 

05 = 3580.98046875 , = -63982.8134766 , 07 = 1329733.72705 , 

as = -31448214.6928 , Og = 833541603.263 , aio = -24478940702.8 . 
Since here Oq = 1/2, the reduced function f{g), for which /(O) = 1, has to be defined as 
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with the corresponding change of the expansion coefficients in series (2). The strong-couphng 
hmit of the ground state energy is 

E{g) ~ 0.667986 g^^^ {g ^ oo) . (20) 
For the self-similar factor approximant E2{g), we have 

Ai = 17.5973 , A2 = 5.3122 ; 

m = 0.0221 , n2 = 0.2091 . 

The large-g limit is 

E;{g) ~ 0.7554 g^-^^^^ {g ^ 00) , 

Comparing the latter with the asymptotic form (20), we see that the error of the amplitude 
here is 13% and that of predicting the power is —31%. 
In the case of E^{g), we get 

Ai = 26.7402 , A2 = 12.4688 , A3 = 3.8380 ; 

ni = 1.8017 X 10-^ , 712 = 0.0547 , rig = 0.2005 . 
The strong- coupling limit gives 

£;*(^)c. 0.7562 /-^s™ (^^00). 

Here the error of the amplitude is again 13%, while that of the power is —23%. If we compare 
E^{g) with the Pade approximant P[3/3](g'), providing the best accuracy at finite g, then we 
see that the latter is wrong at large g, tending to a constant, as g ^ 00, instead of increasing. 
For El{g), we find 

A^ = 36.27784 , A2 = 20.11620 , A^ = 9.69400 , A4 = 2.98137 ; 

ni = 1.05387 x 10"^ , = 8.41555 x 10"^ , ng = 0.08068 , = 0.18273 . 

In the strong-coupling limit, we have 

£;:(^)c^ 0.75229^™ (^-00). 

The errors are around 13% for the amplitude and —18% for the power. 
The approximant E^{g) contains 

Ai = 46.160836 , A2 = 28.182719 , A3 = 16.332980 , 

A4 = 7.924575 , A5 = 2.422310; 
ni = 4.992768 x lO"'' , = 8.787692 x 10^^ , = 0.018549 , 

714 = 0.097564 , rig = 0.164673 . 



15 



The strong-coupling limit yields 

Now the error of the amphtude is 12% and that of the power is —15%. There is a monotonic 
numerical convergence to exact results. 

The accuracy of the approximants El{g) can be essentially improved with the help of the 
crossover condition (7). Invoking the latter, we get for E^i^g) 

= 9.457715 , A2 = 0.869296 ; 

ni = 0.140915 , 712 = -0.192419 ; 

- A^-" = 0.667986 , ni + ria = ^ . 

For E^{g), we find 

Ai = 18.112324 , A2 = 5.896948 ; A3 = 0.483681 ; 
ni = 0.019178 , n2 = 0.184859 , = 0.129296 ; 

^ A1' Al^ Al^ = 0.667986 , m + 712 + rig = - . 
In the approximant El{g), the controllers are 

Ai = 27.098877 , A2 = 12.848678 ; A3 = 4.253915 , A4 = 0.322173 ; 
ni = 1.616412 X 10~^ , 713 = 0.049729 , rig = 0.184752 , m = 0.097236 , 

with condition (7) being valid. 
And for El{g), we obtain 

Ai = 36.555891 , A2 = 20.393831 , A3 = 9.992627 , 

Ai = 3.302859 , A5 = 0.236809 ; 

ni = 9.656933 x 10"^ , 712 = 7.776204 x 10"^ , rig = 0.075166 , 

714 = 0.172047 , ris = 0.078248 , 

again with condition (7) exactly satisfied. 

6 Critical points and critical indices 

The self-similar factor approximants (5) look very appropriate for describing critical phe- 
nomena, where various thermodynamic functions usually possess power-law behaviour. The 
direct application of Fade approximants [1] or their variants [33-35] does not allow for a 
correct description of critical indices. But the factor forms, analogous to approximants (5), 
should provide an accurate account of typical thermodynamic characteristics in the vicin- 
ity of critical points [36]. Note that critical phenomena exist not only in the models of 
statistical mechanics and condensed matter physics but also in the models of rupture and 
materials failure [37-40], earthquakes [41,42], and stock market crashes [43-45]; all such crit- 
ical phenomena having many common features [46,47]. Below we show that the self-similar 
factor approximants present, probably, the best possibility for characterizing different critical 
phenomena. 
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6.1 Logarithmic divergence in elliptic integral 

This example illustrates that even when the divergence in a critical point is not of power 
law but logarithmic, it can, anyway, be reasonably well described by approximants (5). For 
this purpose, we consider the elliptic integral 



TT Jo 



2 r/'^ dx 



f211 



. sm X 



which diverges logarithmically at tc = 1- The logarithmic divergence implies that if one 
defines the critical index as 

a^^hmji-g| ln/(i), (22) 

then it has to be close to zero. Expanding integral (21) in powers oft, one gets the asymptotic 
series of type (2), with the coefficients 

1 9 25 1225 



4' ^64' " 256' " 16384' 
3969 53361 184041 



65536 ' 1048576 ' 4194304 ' 

41409225 _ 147744025 _ 21334223721 

~ 1073741824 ' " 4294967296 ' ~ 68719476736 ' 

For the approximant /2(t), we have 

= -0.9640 , A2 = -0.3485 ; 

ni = -0.2218 , 712 = -0.1038 . 

Divergence happens at a critical point tc, with a critical index, defined in Eq. (22), which 
are 

= = 1.0373 , q; = Inil = 0.222 . 
\Ai\ ' ' ' 

The critical point is defined with an error 3.7%. 
For f*{t), we find 

Ai = -0.9839 , A2 = -0.6311 ; A3 = -0.1767 ; 

ni = -0.1938 , 712 = -0.0721 , 713 = -0.0783 ; 
Divergence is characterized by 

tc=^ = 1.0163, a = |7ii| = 0.194 . 

\Ai\ ' ' ' 

Here, the error of tc is 1.6%. 
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The controllers of the approximants fl{t) are 

Ai = -0.99112 , A2 = -0.77048 , A3 = -0.41456 , 

ni = -0.17705 , 712 = -0.05864 , = -0.05405 , 

The critical point and index are 

1 



A4 = -0.10519 ; 
714 = -0.06587 . 



= 1.0090 , 



a — 72,1 



0.177 



The error of tc is 0.9%. 
For f^{t), we obtain 

Ai = -0.9944 , A2 = -0.8449 , A3 = -0.5782 , 

A4 = -0.2867 , A5 = -0.0692 ; 
ni = -0.1657 , 7i2 = -0.0508 , 713 = -0.0438 , 
714 = -0.0452 , 7i5 = -0.0582 . 
The critical characteristics are 

1 



1.0056 , 



a — \ni\ — 0.166 . 



Now, the error of tc is just 0.6%. Both the critical point and index show good numerical 

convergence to exact results. The best Fade approximant, based on the same number of 
perturbative terms, is much less accurate giving tc = 1.117 and a = 1; the error of tc being 
12%. 



6.2 Critical indices from Wilson expansion 

The scalar quantum field theory with 99^-interaction is an ideal ground for testing new 
methods in perturbation theory [2,48]. Much efforts has been invested into studying critical 
phenomena by means of the 0(7i)-symmetric (/j^-theory. A very popular problem consists in 
calculating critical exponents by field theoretic rcnormalization group techniques in d — A — e 
dimensions, resulting in the Wilson e-expansions [2,49-51]. To illustrate a good accuracy of 
self-similar factor approximants in summing such e-expansions, we consider here the critical 
index 1/ for the correlation length ^{T), defined as 

u=- lim (T-Tc)-^ lnC(T) , (23) 

where T is temperature. At the critical temperature Tc, the correlation length diverges as 
^ ~ (T — Tc)'". We shall use the expansion [49] for the function 
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Prom here, one can get the critical index in three dimensions 



"=27(1) ^'-'^ 

or in two dimensions 

Function (24), presented as expansion (2) in powers of e, in the single-component case 
possesses the coefficients 

oi = -0.1665 , 02 = -0.05865 , ag = 0.06225 , 04 = -0.1535 , 05 = 0.4755 . 

The straightforward usage of such expansions does not give reliable results, since the £-series 
are asymptotic and divergent. For instance, the critical index 

^fe = ^7^ (£^4-d), (27) 

defined through a /c-order perturbative polynomial fk{s), as in Eq. (2), for d = 3 becomes 

i/i = 0.6 , U2 = 0.645 , 1/3 = 0.597 , = 0.731 , = 0.431 (d = 3) . 

The sequence {i^k}, evidently, does not converge. 
For the factor approximant /I (£), we obtain 



Ai = 3.018517 , A2 = -0.439949 ; 

I] 

The corresponding index 



ni = 6.874957 x 10"^ , 712 = 0.425622 



in = 0.634 for = 3 and = 1.216 for d = 2. 

We know that for the two-dimensional Ising model, when e — 2, the index u — 1. 
Imposing this restriction, we find for d = 3 the index 1/3 = 0.629. This is to be compared 
with the Borel-transform based calculations for the e-expansion, yielding 1/ — 0.628 ± 0.001 
[52] and p = 0.629 ± 0.003 [53] and also with the lattice numerical results u = 0.631 ± 0.002 
[2]. 

6.3 Critical characteristics from high-temperature expansions 

High-temperature scries expansions are often used in studying critical phenomena. Here, we 
shall consider the so-called (2 + l)-dimensional Ising model, whose Hamiltonian reads 

^ = E (1 - - E - ^ E '^f ' (29) 

i <ij> i 
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in which erf are PauU matrices; couphng parameter; /i, magnetic field; summation is over 
a 2-dimensional spatial lattice; < ij > denotes nearest-neighbour pairs of sites. This model 
is in the same universality class as the three-dimensional Ising model. High-temperature 
expansions for this model are equivalent to the weak-coupling expansions in powers of g. 
Consider, first, the mass gap 

F{g) = E^{g) - Eo{g) , (30) 

being the energy difference between the first excited and ground state levels, at zero magnetic 
field. For the reduced quantity 

fig) ^ ^ , m = 2 , (31) 

the series expansions in powers of g are of type (2). 

The coefficients in series (2) for a square lattice are [54] 

3 9 11 

ai — — 2 , a2 — — 1 , Os — — - , 04 — — - , — — — , 

ae = -10.253906 , 07 = -28.849609 , Og = -57.418152 , 

Qg = -175.053360 , aw = -365.267989 . 

When g approaches the critical value gc, the mass gap tends to zero, as F ~ {gc — gY, with 
the critical index 

u=- lim (g^-g)-^ ln/(^) . (32) 

Note that the coefficients in series (2), for the present case, are of the same sign, which 
is known to be a difficult situation for summation. The absolute values of these coefficients 
are quickly increasing, which manifests the series divergence. 
For the self-similar factor approximant /Kg'), we find 

Ai = 4.03553 , A2 = -3.03553 ; 

ni = -2.4905 x 10"^ , ^2 = 0.65555 . 
The critical coupling and critical index are 

gc = = 0.32943 , i/ = 712 = 0.656 . 

For the factor approximant fs{g), we have 

^1 = -3.03055, ^2 = 2.34432 + 1.17108 i; A3 = A; ; 

ni = 0.65826 , = (1.57037 + 5.33043 i) x 10"^ , ris = , 

which gives the critical characteristics 

gc = ^ = 0.32997 , u ^ m ^ 0.658 . 
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Approximant fl{g) possesses 

Ai = -3.03764 , A2 = 2.80311 , A3 = -0.85714 + 1.50878 i, = A* ; 

ni = 0.6514 , 712 = -7.37667 x 10"^ , ns = (-7.68487 + 4.56213 i) x 10"^ , 714 = . 
This results in 

Oc = —-- = 0.32920 , = ni = 0.651 . 

\Ai\ 

Finally, for f^{g), we obtain 

Ai = -3.04301 , A2 = 2.85769 , A3 = -2.08221 , 

A4 = 0.05019 + 1.48535 i, A5 = Al ; 
ni = 0.64147 , ^2 = -6.64215 x lO'^ , = 0.02454 , 

714 = (-1.27747 + 7.47775 i) x lO'^ , 715 = 7i^ . 
The corresponding critical characteristics are 

gc = ^ = 0.32862 , u = m = 0.641 . 

The found results can be compared with other calculations, summarized in Ref. [54], where 
the critical coupling is in the interval 0.328 < gc< 0.329, while the critical index lies in the 
range 0.629 <v< 0.646. This is in good agreement with our estimates. 
Let us also consider the susceptibility 

1 d'^Eo , , 

^^^^ = " TV ^ ^=0 ' ^^^^ 

which behaves as xis) ~ {gc~g) "' in the vicinity of the critical coupling g^, with the critical 
index ^ 

7= Urn (gc-g)— In xig). (34) 

g^Qc-o dg 

Expanding the susceptibility (33) in powers of g, one gets the series of type (2). For the 
triangular lattice the coefficients are 

Qi = 6 , 02 = 32.95 , 03 = 166.5 , 04 = 843.046875 , 05 = 4218.416666 , 

ag = 20941.023004 , aj = 103361.512587 , = 507986.371687 , 

ag = 2488222.50870 , = 12155136.2137 . 

Again, all coefficients are of the same sign and fastly increasing in their values. Such a series 
is divergent for any finite g. 

For the factor approximant X2(5')) have 

Ai = 6.22771 , A2 = -4.58786 ; 
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m = -0.03523 , 712 = -1.35562 . 
The critical coupling and index are 

9c= 0.21797 , 7 = |n2| = 1.356 . 

In the next order, we get xliQ)^ with the controllers 

Ai = -4.81702 , A2 = 1.6542 + 2.20379 i ; A3 ^ ; 
ni = -1.18552, 712 = 0.22051 + 0.09988 i , 713 = . 

Prom here, 

= ^ = 0.20760 , 7= |7ii| = 1.185 . 
1^1 1 

For xlid)^ we find 

Ai = -4.76379 , A2 = -0.89137 , A3 = 0.79806 + 1.56616 i , A4 = A^ ; 

7ii = -1.25533 , 7i2 = 2.32460 , 713 = 1.06769 - 0.12380 i, 714 = 7^; , 
which yields 

5c = ^ = 0.20992 , 7 = \ni\ = 1.255 . 
I^il 

And the approximant xtid) contains 

Ai = -4.76723 , A2 = -2.66216 , A3 = -1.30481 , 

A4 = 0.84329 + 1.55612 i, A5 ^ AI ; 
m = -1.24715 , 7i2 = -0.10924 , 713 = 1.54693 , 

714 = 1.02499 - 0.01717 i , n^ = nl. 
The corresponding critical characteristics are 

5e = ^ = 0.20977 , 7= |7ii| = 1.247. (35) 
l^il 

These results are in good agreement with numerical calculations, summarized in Ref. [54], 
according to which the critical coupling for the triangular lattice is in the interval 0.20972 < 
Qc < 0.20980 and the critical index is in the range 1.236 < 7 < 1.250. 
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7 General discussion and possible variants 



The self-similar factor approximants, as is illustrated above by a variety of examples, provide 
a powerful tool for defining effective sums of divergent power series. The latter in the 
majority of applications are asymptotic. However, the presented method is not limited to 

only asymptotic scries. The sole thing that is required is that the pcrturbative form (2) be 
a finite power series. If all coefficients a„ in Eq. (2) do not depend on the approximation 
order A; = 1, 2, . . ., then the difference 

^fk{x) = fk{x) - fk-i{x) 

composes as asymptotic sequence {Afk{x)}, where Afk{x) ~ x^. But we also may 
consider the case [5], when an iterative procedure results in the approximations 

k 

fk{x) = 1 + ^ a„fe , 

n=l 

with different coefficients a„fc 7^ dnk-i- In such a situation, the difference Afk{x) does not 
constitute an asymptotic sequence {Afk{x)}. Nevertheless, the method of self-similar factor 
approximants is applicable to this case as well. 

The mathematical structure of formula (5) resembles, in some sense, such geometric 
structures as fractals [55,56]. With regard to functions, one usually says that a function 
possesses fractal properties if it manifests a kind of scaling. In particular, it can be an 
asymptotic scaling [56,57], arising in a limit of the function variable. In that sense, formula 
(5) displays the property of asymptotic scaling 

hm = A-^S = E , (36) 

where A is a scaling parameter. We also may notice another scaling property of the factor 
form (5). To this end, introducing functions 

ip,{x) = {l + A, x)"^ , (37) 

we can present Eq. (5) as 

k 

f:ix)^l[^,{x). (38) 

i=l 

Then, the property of functional scaling 

fki^) ^ (n fki^) ' ^ii^) ^ ^ii^) ' (39) 

is valid. 

It is worth emphasizing that the factor approximant (5) is principally different from 
the continued- function representation mentioned by Bender and Orszag [58]. According to 
their scheme, approximants (3k{x) for a function f{x) could be constructed in the following 
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way. Let us take an arbitrary function P{c,x), where c is a parameter. Then a /c-order 
approximation to the sought function f{x) is given by 

Pk{x) = P (Ci, XP (C2, . . . , xP{Ck, X)) ...) . 

However, such a manner of acting has several essential drawbacks. First of all, this scheme 
is too much ambiguous, since there exists an infinity of functions that could be accepted as 
f3{c, x). Contrary to this, the form of the factor approximants (5) is not arbitrary postulated 
but is derived from the given perturbative scries by means of the self-similar approximation 
theory. Second, Eq. (5) does not have the form of a continued function Pk{x). Really, even 
if we suppose that /3(c, x) = {1 + cx)", then 

Pkix) = (1 + CiX (1 + C2X{1 + ... + Ck X)")" ...)", 

which has nothing to do with formula (5). The latter is not a function of the same func- 
tion but a product (38) of different functions (37). This product is not only different from 
the above /3k{x) but also is more general. Last but not least, the sequence {(3k{x)} of 
the continued-function representation is, generally, speaking, divergent. This can be easily 
demonstrated for the above mentioned case of (3{c,x) resulting in the power- law representa- 
tion Pk{x)- The latter, at large x, arc of the type 




{x — > oo) , 



which shows that the subsequent approximations l3k{x) differ significantly from each other. 
Thus, f3i{x) ~ x", f32{x) ~ a;^", and so on, each f3k{x) ~ x"^ having principally different 
dependence on x. It is evident that such a sequence cannot converge. 

The self-similar factor approximant (5) of order k requires the knowledge of 2k nontrivial 
(except the first term equal to one) perturbative terms of series (2). Thus, for constructing 
fk{x), we need an even number of terms. A natural question would be: how could we use 
an odd number of such terms? There are two straightforward ways of doing this, which can 
be called additive and multiplicative. 

The additive way is as follows. The perturbative form (2) can always be rewritten as 

k-l 

fkix)^l + aixl[{l + Cix) . (40) 

i=l 

By applying the self-similar renormalization (4) to the factors in the product of Eq. (40), 
we come to the approximant 

r,^{x) = l + a,xll{l + B,xr' , (41) 

i=l 

containing the additive unitary term. The approximant (41) requires the availability of an 
odd number of terms in series (2). If the large-a; behaviour (6) is known, it can be employed 
for approximant (41) as the crossover condition 

k k 

ail[Br=foo, J2^i + 1^(3. (42) 
i=i i=i 
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We have analysed the usage of approximants (41) to the examples considered above. Some- 
times these approximants fk+{x) are not bad, being of comparable accuracy with /^(a;) or 
even better. This, e.g., is the case for defining the critical index (23) from the Wilson e- 
expansion, for which we find ^'2_|_ = 0.628 for d = 3 and i>2+ — 0.914 for d = 2. However, 
for the majority of cases, the accuracy of fk+{x) is essentially worse than that of fl{x). On 
occasion, this accuracy can be improved by considering the inverse function f^^{x), instead 
of f{x). This, for instance, gives good results for the structure factors of Section 4. The 
crossover condition (42), if available, also may certainly cure the behaviour of f^^{x). For 
example, this yields good approximations for the ground-state energy of the anharmonic 
oscillator of Section 5.2. But, in general, the approximants fl_^{x) are inferior to fl{x). 

As an alternative, we could try the multiplicative way of extracting information from odd 
numbers of perturbative terms. For this purpose, we may present Eq. (2) as 

/fe(a;) = (l + ^x)n(l + 6ix). (43) 
1=1 

Accomplishing the self-similar renormalization of the right-hand-side factors in Eq. (43), 
except the first factor (1 -|- Ax), we define the self-similar approximant 

/;+i/2(^) = (1 + Ax) n(l + A x)-^ (44) 

1=1 

In the starting trivial case, 

/i72(^) = /o(^) = l + «i^- (45) 
The following nontrivial approximants are 

n/2{x) = {l + Ax){l + A,x)''\ 

f;/,{x) = {l + Ax)il + A^x)^^{l + A,x)^' , (46) 

and so on. Each fk+i/2{^) requires for its definition 2A: + 1 terms of series (2). 

The main difference in employing approximants fk+i/2{^)j ^ compared to fk{x), is in 
the following. When determining the related control parameters by means of the order- 
through-accuracy matching, for the former approximants we may get nonunique solutions, 
while for the latter the solutions are unique. Different solutions for the control parameters 
yield several possible variants of fk+i/2{^)- '^'c have checked, all of these several variants, 
generally, have reasonable accuracy, being close to each other. However, dealing with several 
approximants of the same order complicates the calculational procedure. In this way, the 
self-similar factor approximants /^(x) look to be more preferable, compared to both f^^{x) 
and fl+i,2{x). 

In all cases, we have considered, the accuracy of fil{x) improves with the increasing order 
k. This implies the existence of the so-called numerical convergence for the sequence {/^ (a;)}. 
It is, of course, impossible to guarantee that the latter sequence always converges to the exact 
function f{x). This kind of uncertainty, is actually, common for practically any sequence of 
numerical approximants. To render the procedure uniquely defined requires imposing on the 
calculational algorithm additional conditions making the latter single- valued [59]. 
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More information of the self-similar factor approximants fk{x) can be extracted by in- 
voking the probabilistic analysis [60]. To this end, we introduce the multiplier 



Mk{x) 



(47) 



Taking account Eq. (8), we get 



1 d 
ai dx 




(48) 



The probabihty of fl{x) is defined as 



M^\x)\ 



(49) 



Ef=i \M-\x)\ ' 



where is the total number of available approximants f^ix), with k = 1,2, . . . , N. Having 
defined these quantities, we now may find the most probable approximants, corresponding 
to the maximal probability (49). Also, we can calculate a weighted approximant 



This averaging is especially useful when the subsequent approximants fj{x) and f*_^_i{x) 
display large variations. The latter become smoothed by the averaging procedure (50). For 
example, for the critical index of the correlation length, discussed in section 6.2, defined by 
this averaging procedure, we find u = 0.630. 

Additional possibilities arise, when, for the same problem, there exist expansions for 
several different functions. Say, for simplicity, there are two functions f{x) and g{x), for 
which the expansions fk{x) and gk{x) are available at x — > 0. Then, in addition to the factor 
approximants /^(x) and glix), we may, considering some combinations Fk{fk{x),gk{x)) of 
the series fk{x) and gk{x), construct the approximants F^{x). For instance, we can consider 
the sums fk{x)+gkix), products fk{x)gk{x), or the ratios fk{x)/gk{x), or other combinations 
of the given series. A more detailed consideration of such possibilities will be done in future. 

Finally, we would like to emphasize that the accuracy of the self-similar factor approx- 
imants (5), for the examples considered in the present paper, is always essentially better 
than that of the best Fade approximants, constructed with the same number of perturbative 
terms of series (2). This fact is easy to understand, since the factor approximants (5) have 
a much more general form as compared to Fade approximants, including the latter just as 
a narrow particular class. As is explained in Section 2, the self-similar factor approximants 
can represent rational, irrational, and transcendental functions. 
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Figure Captions 



Fig. 1. Difference of the self-similar factor approximant and function (13). 

Fig. 2. Specific heat (14) (solid line), the approximant C|(x) (dashed line), and the 
best Pade approximants P[a/2]{x) (dashed-dotted line) and P[5/i](a;) (dotted line) versus 1/x. 

Fig. 3. Percentage errors of the approximant f^{x) (solid line) and the two-point Pade 
approximant P[3/2](a;) (dashed line), as compared to the exact function (15). 

Fig. 4. Structure factor of branched polymer S{x) (solid line) and its factor approximant 
5*2 (x) (dashed line). 

Fig. 5. Luminiscent intensity f{z) (solid line), approximant f^{z) (dashed line) and the 
best Pade approximant P[3/3](a;) (dashed-dotted line). 
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